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Abstract 

In three-dimensional quantum electrodynamics (QED 3 ) with mas¬ 
sive gauge boson, we investigate the Dyson-Schwinger equation for the 
fermion self-energy in the Landau gauge and find that chiral symmetry 
breaking (CSB) occurs when the gauge boson mass f is smaller than 
a finite critical value £ c „ but is suppressed when £ > f,cv • We further 
show that the critical value £ c „ does not qualitatively change after 
considering higher order corrections from the wave function renormal¬ 
ization and vertex function. Based on the relation between CSB and 
the gauge boson mass £, we give a held theoretical description of the 
competing antiferromagnetic and superconducting orders and, in par¬ 
ticular, the coexistence of these two orders in high temperature super- 
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conductors. When the gauge boson mass £ is generated via instanton 
effect in a compact QED 3 of massless fermions, our result shows that 
CSB coexists with instanton effect in a wide region of £, which can be 
used to study the confinement-deconfinement phase transition. 
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1 INTRODUCTION 


Chiral symmetry breaking (CSB) has been an active research field in particle 
physics for over forty years since Nambn [1] used this idea to generate fermion 
mass in a four-fermion model. One fascinating characteristic of CSB is that 
it can generate fermion mass via the fermion-antifermion condensation me¬ 
diated by strong gauge field without introducing additional Higgs particles 
which till now have not been found. The most conclusive evidence for the 
existence of CSB is provided by the phenomenology of strong interaction, 
and CSB is widely believed to account for the pions. However, despite the 
vast amount of theoretical work on CSB, till now it is not clear whether it 
can be derived from quantum chromodynamics (QCD), primarily due to the 
complex structure of the SU(3) gauge field. To gain valuable insight of CSB 
before we can treat it completely, it is very suggestive to study some model 
that is similar to QCD while is simpler than it. Three-dimensional quantum 
electrodynamics (QED3) is just such a model, and it has attracted intense in¬ 
vestigations [2-9] in the past twenty years. QED3 was shown to exhibit CSB 
[5-9] and confinement [9], while at the same it is simple enough to be treated 
with high accuracy. Besides, it has been used to model the physics of many 
planar condensed matter systems such as high temperature superconductors 
[10-18] and fractional quantum Hall systems [19]. 

The breakthrough in the research of CSB in QED3 was brought by a paper 
of Appelquist et al. [5] who found that CSB occurs when the flavor of massless 
fermions is less than a critical number N c . They arrived at this conclusion by 
analytically and numerically solving the Dyson-Schwinger (DS) equation of 
the fermion self-energy to the lowest order of 1 /N expansion. Later, extensive 
analytical and numerical investigations [6-8] showed that the nature of CSB 
in QED 3 remain the same after including higher order corrections to the DS 
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equation. 

The above result holds when the gauge boson is massless and but is 
expected to change when the gauge boson has a finite mass. CSB is a low- 
energy phenomenon because (2+l)-dimensional U(l) gauge field theory is 
asymptotically free [20] and only in the infrared region the gauge interaction 
is strong enough to cause fermion condensation. This requires the fermions 
be apart from each other. However, when the gauge boson has a finite mass 
it can not mediate a long-range interaction. Intuitively, a finite gauge boson 
mass is repulsive to CSB which is achieved by the formation of fermion-anti¬ 
fermion pairs. Thus it is very interesting to study whether CSB can occur in 
the presence of a finite gauge boson mass. 

CSB is believed to be a nonperturbative phenomenon and hence calcu¬ 
lations based on perturbative expansions are incapable of establishing its 
existence. We will study CSB by means of solving the nonlinear DS equation 
for the fermion self-energy. Assuming that A(p 2 ) = 1 based on naive 1/N 
expansion, we get a single integral equation of the gap function E(p 2 ). CSB 
is signalled by the appearance of a squarely integrable nontrivial solution. To 
solve the DS equation, we will use bifurcation theory and parameter imbed¬ 
ding method, which not only avoids the convergency problem that usually 
appears in iteration method but also can help us distinguish the different 
bifurcations points. After solving the DS equation, we find that the massless 
fermions can acquire a finite dynamically generated mass when the gauge bo¬ 
son mass £ is smaller than a critical value £ c . To testify the robustness of our 
result against the effect of A(p 2 ), we will work in a non-local gauge in which 
A{p 2 ) = 1 and the vertex function can be replaced by the gamma matrices 
safely. We will show that when the wave function renormalization A(p 2 ) 
is included, the result we derive in the Landau gauge remains qualitatively 
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unchanged. 

Our study on the influence of gauge boson mass on the fate of CSB not 
only is of theoretical interests but can be used to understand important physi¬ 
cal phenomena. Actually, starting from the concept of spin-charge separation 
proposed by Anderson [21], the effective low energy theory of high temper¬ 
ature superconductors is a U(l) gauge theory [10-13]. Superconductivity is 
achieved when the charge carrying holons Bose condensate into a macro¬ 
scopic quantum state, which generate a finite mass to the gauge boson via 
Anderson-Higgs mechanism. The low energy spin fluctuations are captured 
by the two-component fermions, which are originally massless since they are 
excited from the d -wave gap nodes [11]. On the other hand, CSB is known 
to correspond to the long-range antiferromagnetic (AF) order, which can be 
seen from the behavior of AF spin correlation function at low momentum. 
If we use CSB to describe the AF order and use the gauge boson mass to 
describe the superconducting (SC) order, our result then leads to a competi¬ 
tion between the long-range AF order and the long-range SC order, which is 
one of the most fundamental issues in modern condensed matter physics. As 
a compromise of this competition, when the mass of gauge boson is less than 
its critical value £ c but is finite there is a coexistence of these two orders in 
the bulk superconductors. 

If the U(l) gauge field is compact in the meaning that the vector poten¬ 
tial has a periodicity, then it acquires a finite mass via the instanton effect. 
Furthermore, permanent confinement of static charges is present when the 
instanton effect is important. The influence of additional matter fields, es¬ 
pecially massless fermions, on the permanent confinement is an unsolved 
problem. The relation between CSB and the gauge boson mass obtained 
in this paper is very helpful in studying the confinement to deconfinement 
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phase transition driven by the couling of massless fermions. 

The physical applications of CSB in the presence of a gauge boson mass to 
high temperature superconductors have been reported in a Letter [16]. In this 
paper, we provide the related held theoretical technique in details. In Sec. II, 
we derive the DS equation in the presence of a finite mass of the gauge boson 
in the Landau gauge. We then choose to solve the nonlinear DS equation 
by means of bifurcation theory and parameter imbedding method. Sec. Ill 
is devoted to the elementary knowledge of bifurcation theory that will be 
used in this paper and the detailed calculation steps of parameter imbedding 
method. In Sec. IV, we consider the higher order corrections to the wave 
function renormalization and show that these corrections do not change our 
result noticeably. In Sec.V, we give a thorough discussion of the competing 
orders in high temperature superconductors from a held theoretical point of 
view. In particular, we emphasize the necessity for nonperturbative effect in 
getting an AF spin correlation that is consistent with experiments. In Sec. VI, 
we discuss the instanton effect on CSB in compact QED 3 . The calculation 
of AF spin correlation function in the CSB phase is given in the Appendix. 

2 DYSON-SCHWINGER EQUATION IN THE 
LANDAU GAUGE 

The three-dimentional U(l) gauge theory of massless fermions is 

1 N — 

£ = -A + Y ( d H - ia n) 7mU, (1) 

( 7=1 

where the ferrni held ^ is a 4 x 1 spinor. The 4 x 47 ^ matrices obey the 
algebra, { 7 ^, 7 *,} = 2 8^. 

The full fermion propagator is 

G-\p) = i'y ■ pA (p 2 ) + E (p 2 ) , (2) 
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where A(p 2 ) is the wave-function renormalization and E(p 2 ) the fermion self¬ 
energy. The DS equation for the full fermion propagator in momentum space 
is given by 

G~\p) = Go\p) ~ I ■^ E 'y fi G(k)r u (p, k)D^(p - k), (3) 

where T u (p, k ) is the full vertex function and D^ u (p — k ) is the full photon 
propagator. Gq 1 {p) is the bare propagator of the massless fermions. Sub¬ 
stituting the propagator (2) into Eq.(3) and taking trace on both sides, we 
obtain the equation for E(p 2 ) 

s (p 2 ) = ~\J -j7^ Tr b » G (k)T v (p, k)D^(p - k)}. (4) 

Multiplying both sides of Eq.(3) by 7 -p and then taking trace on both sides, 
we obtain the equation for A(p 2 ) 

A(p 2 ) = 1 + 4^2 / ?M(*7 ■ fc)7 y.G{k)V v {p, k)D^(p - k)}. (5) 

If the DS equation for E(p 2 ) has only vanishing solutions, the fermions 
remain massless and the Lagrangian (1) respects the chiral symmetries ^ —» 
exp(?' 073 !5 )^, with 73 and 75 two 4x4 matrices that anticommute with 7 M 
(/j = 0,1,2). If the DS equation for E(p 2 ) develops a squarely integrable 
nontrivial solution [23-25], then the originally massless fermions acquire a 
finite dynamically generated mass which breaks the chiral symmetries. 

We starts from a general gauge boson propagator 


with 


DM = D T (q 2 ) S^-g(q 2 )^ 


D T (q ) = q 1 +7T(g ) 


( 6 ) 

(7) 
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where g(q 2 ) is a gauge-fixing parameter that depends on 3-momentum and 
£ is the mass of the gauge boson. We use q to denote the gauge boson 
momentum, and we have q 2 = (p — k ) 2 = p 2 + k 2 — 2pkcos9. tt( q 2 ) is 
the vacuum polarization of the gauge boson, which was included initially to 
overcome the infrared divergence. If we includ only the one-loop diagrams 
for massless fermions, we can write the vacuum polarization as 

^ ^ m (8) 

Then we have 


Dr\q 2 ) =qMq 2 )+e 

N 

= 7(9 + '')), 


with 


(9) 


( 10 ) 


<1 = 8 C//V, 

since at low momentum 7 r(g 2 ) S> 1 [5]. 

As the lowest-order approximation, we neglect the wave function renor¬ 
malization A(p 2 ) and adopt a massive gauge boson propagator in the Landau 
gauge ( g(q 2 ) = 1 ) as follows 

(P~ k)„(p- k) v \ 


D,,J.P ~ k) = 


8 „ 


( 11 ) 


N(\p — k\ + rj) \ ^ ( p-k) 2 )' 

Further, we use the bare vertex, i.e., T^{p,k) = 7 ^, which is usally called 
quenched planar approximation. Then the DS equation becomes 

d 3 k ^ D ^(jp — k)T 1 (k 2 )Y 


3 (2vr) 3 k 2 + Z 2 (k 2 ) 

Now we can insert the propagator (11) into the DS equation (11), then 


( 12 ) 




dk 


kZ(k 2 


dz- 


k 2 + YA(k 2 ) J-i \p-k\ + r) 


(13) 



where z = cos 9. After performing the integration with respect to z and intro¬ 
ducing an ultraviolet cutoff A we finally arrive at the following DS equation 


£(p 2 ) 


M 

A / dk 
Jo 


kE(k 2 ) 
k 2 + E 2 (A; 2 ) 


1 

x- 

P 



( p+k+p\\ 

\\p-k\ + v)) ’ 


(14) 


where A = 4/Abr 2 serves as an effective coupling constant. 

If we do not introduce an ultraviolet cutoff (A —> oo), the critical behavior 
of Eq.(14) is completely independent of 77 , as can be easily seen by making 
the scale transformation, p —> p/ 77 , k —> k/r 7 and E —> E/ 77 . We can destroy 
this scale invariance by introducing an ultraviolet cutoff A. 

Before we go into the techniques of dealing with the DS equation, we 
would like to discuss one subtle issue. If the DS equation has only trivial 
solutions, the fermions remain massless and the chiral symmetries are not 
broken. However, not all nontrivial solutions lead to CSB. It is well known 
that the breaking of a chiral symmetry is always accompanied by a Goldstone 
boson, which is a pseudoscalar bound state composed of a fermion and an 
antifermion. If CSB happens, there should be a nontrivial solution for the 
Bethe-Salpeter (BS) equation of this bound state. In addition, the bound 
state wave function must satisfy a normalization condition, which can be 
converted to a sufficient and necessary condition [23-25] for the nontrivial 
solutions of the DS equation to signal CSB. It gives a constraint on the form 
of E(p 2 )/A(p 2 ) as follows [25] 



g 2 £V) 

q 2 A 2 (q 2 ) + S 2 (g 2 ) 


finite. 


(15) 


It is easy to see that in order to satisfy this condition E(p 2 )/A(p 2 ) must 
damp more rapidly than p ~ 1//2 in the ultraviolet region (p —> 00 ). The mass 
function obtained by Appelquist et al. [4] satisfies this condition and hence 
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the nontrivial solutions of the DS equation in QED 3 corresponds to true CSB 
solutions [25]. On the other hand, when an ultraviolet cutoff is introduced, 
the solution E(p 2 ) automatically satisfies the squarely integrable condition. 
We should emphasize that although the nontrivial solutions with an ultra¬ 
violet cutoff all satisfy such a condition, only those solutions that satisfy 
this condition in the continuum limit is taken are physically sensible. In the 
case of four-dimensional QED, although the nontrivial solutions with explicit 
ultraviolet cutoff are squarely integrable they do not satisfy the squarely in¬ 
tegrable condition when we take the continuum limit [23,24], Therefore, the 
CSB solutions obtained in quenched planar QED 4 [26] are not physically 
meaningful solutions because they are not squarely integrable in the contin¬ 
uum limit, or in other words the associated bound state wave functions can 
not be normalized. In QED 3 , it was found that T,(p 2 )/A(p 2 ) behaves like p~ 2 
at p —> oo, which surely satisfies the condition. Since the nontrivial solutions 
in QED 3 are true CSB solutions, we can safely introduce an ultraviolet cutoff 
A without bringing unphysical nontrivial solutions. 

Theoretical analysis implies that the critical fermion number N c of Eq. (4) 
should depend on A/ 77 . To determine when CSB occurs, the DS equation 
should be solved implicitly. The DS equation is an nonlinear integral equation 
and hence is very hard to investigate. However, based on general bifurcation 
theory and parameter imbedding method, we can find the critical fermion 
number and the mass function exactly. The detailed programe of calculations 
is the topic of the next Section. In the rest of this Section, we discuss some 
qualitative properties extracted from the DS euqation (14). 

When the gauge boson has a very large mass, for example 77 A, then 
the DS equation becomes 


£b 2 ) 


8 [* dk k 2 ) 

Nir 2 ri Jo k 2 + T 2 {k 2 )' 


(16) 
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From the momentum dependence of the mass function of fermions, we know 
that actually it is a constant in this limit. Therefore, the DS integral equation 
simplifies to a algebraic equation 

£ arctan (^j = ^ (l - vr 2 //) . (17) 

This equation has no solutions, hence a large enough mass of gauge boson 
prevents the occurrence of CSB. We now consider another limit, i.e., when 
the gauge boson mass is very small. In this limit, the last term in the kernel 
of (14) can be dropped safely, leaving a DS equation that is the same as the 
one studied by Appelquist et al. [5]. Then a very small gauge boson mass 
actually does not affect the critical behavior of QED 3 . This phenomenon can 
be understood if it happens that the critical fermion number N c decreases 
when the gauge boson mass £ increases and it finally approaches zero for very 
large £. The main purpose of our work is to use QED 3 to model condensed 
matter systems where the physical fermion number is 2, which comes from 
the two components of the spin. Based on the tendency of the critical fermion 
number in the presence of a finite gauge boson mass, it is quite reasonable to 
hypothesize that there is a critical value for the gauge boson mass £ cv above 
which CSB is inhibited. To make sure that it is actually the case, we should 
solve the DS equation and find the critical coupling constant A c at which the 
DS equation starts to have nontrivial solutions. 

3 SOLVING DS EQUATION USING BIFUR¬ 
CATION THEORY AND PARAMETER 
IMBEDDING METHOD 

The equation (14) is a Hammerstein type nonlinear integral equation. It 
does not satisfy the conditions of the global eigenfunction theory of nonlin- 
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ear functional analysis, so its global solutions can not be obtained directly. 
However, the local bifurcation theory [27] can help us to find its complete 
solutions by first obtaining a local solution near a bifurcation point and then 
extending its region of validity step by step. This program is most easily 
achieved by parameter imbedding method [28-30], which has proved to be a 
powerful method in studying integral equations. 

In order to obtain the bifurcation points we need only to find the eigenval¬ 
ues of the associated Frechet derivative of the nonlinear DS equation [28,29]. 
Those eigenvalues that have odd multiplicity are the bifurcation points. Mak¬ 
ing Frechet derivative of the nonlinear equation (14), we have the following 
linearized equation 

rA/r] 

E(p 2 ) = A/ dkZ(k 2 )K(p,k ) (18) 

Jo 

with the kernel 

«M = y k (p + k-\p-k \-In ) (19) 

where for calculational convenience we made the transformation p —> p/ 77 , 
k —> k/rj and E —> E/ 77 . The smallest eigenvalue A c of this equation is just 
the bifurcation point from which a nontrivial solution of the DS equation 
(14) branches off. The complex kernel K(p, k ) in the linearized equation (18) 
makes it very difficult to find an analytical solution. 

We now would like to use parameter imbedding method [28,29] to solve 
(18) numerically. To do this, we first analytically continue it in the complex 
plane of A, correspondingly E(p 2 ) also becomes a complex function. It can 
be shown that 

J J \K(x, y)\ 2 dxdy < 00 . ( 20 ) 

Here we use x to denote p 2 , and y to denote k 2 . From the Fredholm integral 
equation theory we know that there exists a resolvent function for the kernel 
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K(x,y) 


R F (x,y, A) = D F (x,y,X)/d F (X), 


( 21 ) 


where the functions D F (x,y, X) and d F (X) are analytic with respect to A. If 
d F ( A) 7 ^ 0, we do not have bifurcations points. The values of A at which 
d F { A) = 0 are the bifurcation points. According to the parameter imbedding 
method, the functions D F (x,y, A) and d F ( A) are related by the differential- 
integral equations 

—d F (X) = - / 17 D f (x, y, X)dx, ( 22 ) 

d X Jo 


d 

dX 


D F (x,y, A) 


d F (X) 


D F (x, y, X)—d F (X) + J ^ 


D f (x, z, X )D f (z, y, X)dz 
(23)' 


with the initial conditions 

d F ( 0 ) = 1 , 

Df(x, y,0) = K(x,y). 


(24) 

(25) 


One remarkable advantage of parameter imbedding method is to convert 
the integral equations with variables x and y to a set of equations in the 
variable A. Correspondingly, the boundary conditions in original equations 
are replaced by two initial conditions, which are easier to treat in performing 
numerical calculations. Now, the functions D F (x,y, A) and d F ( A) can be 
readily obtained by integrating numerically with respect to A. 

We now should choose an appropriate contour C in the complex A-plane 
which contains the minimum A on the real axis at which d( A) = 0. The 
number of zero eigenvalues of the linearized equation (18) inside the contour 
C is (i.e., the zeroes of the function d F ( A)) 

N f = - J) — — ————d F (X)dX. (26) 

2mJcd F (X)dX K J y ’ 
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( 27 ) 


We can obtain the eigenvalues by solving the equations 


VA \ = — i -TTTT-Kd F (X)d\ 

2m Jc d F ( A) dX 


with l = 1, ...,N f . For the present purpose, we only need to know the first 
bifurcation point, hence we let N F = 1. 

For A > A c , the DS equation has nontrivial solutions and the massless 
fermions become massive. The ultraviolet cutoff A is provided by the lattice 
constant and hence is kept fixed. We can obtain the relation of N c and 77 by 
calculating the critical coupling A c for different values of A/ 77 . 

Our numerical result is presented in Fig.(l). The critical fermion number 
N c is a monotonously increasing function of A/ 77 . For small A/ 77 , Al c is smaller 
than physical number 2; so CSB does not occur. When A /77 increases, the 
critical number N c increases accordingly and finally becomes larger than 2 
at about A/r] cv = 100. Then we see that there is a critical value of the gauge 
boson mass £ c „, below which a finite mass is generated for massless fermions 
while beyond which CSB is suppressed. 


4 DYSON-SCHWINGER EQUATION WITH 
HIGHER ORDER CORRECTIONS 

In the last two Sections, we have investigated the DS equation in the Lan¬ 
dau gauge after assuming that A(p 2 ) = 1 to simplify calculations. Although 
this assumption is qualitatively correct, higher order corrections from the 
wave function renormalization A{p 2 ) will alter the critical fermion number 
N c quantitatively. However, including A(p 2 ) makes the DS equations very 
complicated and we should solve consistently two couples of nonlinear inte¬ 
gral equations. Furthermore, according to the Ward-Takahashi identity, we 
can not choose 7 ^ as the vertex function in the presence of wave function 
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renormalization A{jp 2 ). At present, there is no theoretical guidance in deter¬ 
mining the vertex function T M (p, k), and hence one can not give a guarantee 
of the legitimacy of a specific choice of vertex function. Here, to simply 
calculations and partly overcome the embarrassment in choosing the vertex 
function, we introduce a so-called nonlocal gauge [31-34] in which the wave 
function renormalization A{p 2 ) = 1 and the vertex function can be chosen as 

r,(p,k) = ^f(p 2 ,k 2 ) (28) 

with / a function of the fermion momentum p 2 , k 2 . The nonlocal gauge is 
obtained by solving a differential equation. In this gauge, we need only to 
investigate a single equation for E(p 2 ) in studying the chiral phase transition. 

Let us go back to the general massive gauge boson propagator (6). If 
we consider quenched planar approximation of QED 3 , i.e., taking Tl(p 2 ) = 0, 
then the wave function renormalization A{p 2 ) = 1 in the Landau gauge. 
This result is well-known to be exact in QED of dimensions higher than 2. 
In the case of QED 3 , the one-loop vacuum polarization is usually introduced 
explicitly to overcome the severe infrared divergence. In the presence of 
H(p 2 ), wave function renormalization A(p 2 ) does not equal to identity. It 
should be obtained by solving two consistent integral equations of A(p 2 ) and 
£(p 2 ). However, taking advantage of the gauge degrees of freedom of the 
system, we can simplify the DS equations by choosing an appropriate gauge. 
In particular, if we could obtain a gauge parameter g(q 2 ) that satisfies the 
following equation [34] 

9(q2)=DT{z)zdz ~ h (29) 

then we find a gauge in which A(p 2 ) = 1. Further, according to the Ward- 
Takahashi (WT) identity, the vertex function can be chosen as [34] 

r M (p,&)=7/^ ( 30 ) 
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Now the formidable task to solve a pair of integral equations for the wave 
function renormalization A(p 2 ) and the mass function £(p 2 ) is simplified to 
solve a single equation of £(p 2 ) 


£(p 2 ) = 


d 3 k Y>(k 2 ) 


3 — g(q 2 ) D T (q 2 ). 


(27r) 3 k 2 + S 2 (/c 2 ) L 
From D T (q 2 ) and Eq.(29), the integral upon z can be calculated 


r<i 2 8 f 

L ° T ^ zdz = nL 


8 n i 

T -(jaz 


0 z 2+7/ 

16 /l 3 1 2 2 3 (q +7?Y 

= V(3 9 - 2 W + J, 

Then we obtain a nonlocal gauge parameter 

s(« 2 ) = 4(9 + ’ll (l,'f - 7 /r + - <f In (l + -Yl - 1. 


(31) 


r 


n, 


(32) 


(33) 


Substituting this g(q 2 ) into the DS equation (31), after angular integration 
we have 


W) = 


Nti 2 p Jo 

rp+k 

X / dq 
J\p-k\ 


dk 


kT,(k 2 ) 
k 2 + T?{k 2 ) 


2 

3 


V V_ 

q + i7 2 q 


2 S 

V . T 

— + —7 In 

q z q 6 


In deriving this result, we have used the following formula 



f-K 1 rp+k 

/ d0sm0f(q 2 ) = — qdqf(q 2 ). 
Jo pk J\p—k\ 

After integrating (34), we have 

£(/) = A /; dk JfAA 2 -K(p, k, „), 

with 

K (p+ n) = |l> + fc-|p-fc|) ~ | In ( ^ fc,+T? ) 


(34) 


(35) 


( 36 ) 
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T] 

"T 


V 


2 \p — k\ 


p + k 

3 —lnfl + 


\p-k\J 
\p-k r 


T] 


In 1 + 


p + k' 


37) 


7] J 2 (p + k) 2 

Here, A = A/Nn 2 is the effective coupling constant. At the first glance, both 
the third and fourth terms of K(p , k, rf) have singular behaviors like 1/ \p — k\ 
which would cause divergence if k approaches p. However, when \p — k\ —> 0, 
we can make the expansion 


V 


f \p — A;|\ rf f\p — k\ (p — k) 2 


2 |p — k\ 


■In 1 + 


V 


V 


2 if 


2 | p — k | 


+ 


o(|p-fc| 3 ) 


(38) 


Thus the singular terms are exactly cancelled. The same step can be used to 
show that the singular term 1 /(p+q) can also be cancelled exactly. Therefore, 
the kernel K(p, k,r /) is a smooth function on the whole integration region. 

Making Frechet derivative of the nonlinear equation (35), we obtain the 
linearized equation 


E(p 2 ) = A r dkTj(k 2 ) — K(p, k, p) 
Jo pk 


(39) 


with 


K{p, k, rf) = ^(p + k - |p - k\) - ^ In ^ ^ k^+l j 

1 f 1 l \ 

P 2 \p + k \p — k\ J 

+ yyby 111111 + lp - p|) “ WPW (1+p+ km) 

where for calculational convenience we made the transformation p p/p, 
k —> k/p and E —► E/ 77 . 
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Using the steps we presented in the last Section, we can solve the lin¬ 
earized equation (39) to obtain the relation between the critical fermion 
number N c and the mass £ of the gauge boson mass. The numerical result is 
presented in Fig. (2), from which we know that the critical value of the gauge 
boson mass is about A/rj cv = 3.3. Although there is a significant change in 
the critical value £ c „, the result we obtained in the Landau gauge remains 
qualitatively correct. 

5 COMPETING ORDERS IN HIGH TEM¬ 
PERATURE SUPERCONDUCTORS 

Understanding the competing orders in high temperature cuprate supercon¬ 
ductors is one of the most important issues in condensed matter physics. In 
the presence of competing orders, one order parameter prevails when other 
orders are suppressed by some external variables. At half-filling, the cuprate 
superconductor is a Mott insulator with long-range antiferromagnetic (AF) 
order. When holes are dopped into the Cu-0 planes, the material becomes 
a superconductor at low temperatures and the long-range AF order disap¬ 
pears. Hence there is a competition between the AF order and the SC order, 
and as a result of this competition the AF order dominates at zero and low 
doping while the SC order dominates at higher doping. However, even at 
higher doping the AF order also has a chance to appear locally where the 
superconductivity is suppressed by strong external magnetic fields. Recently, 
elaborate neutron scattering [35] and scanning tunnelling microscopy (STM) 
[36] experiments found that the AF correlation is significantly enhanced in 
regions surrounding the vortex cores. In this paper, we will use spin-charge 
separation and CSB to understand the competing orders. 

It has been shown that [10-13] Lagrangian (1) is the effective low energy 
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theory of undoped cuprates which has only fermionic excitations because of 
the presence of a large charge gap. In underdoped cuprates, the electrons 
fractionalize into spin carrying spinons and charge carrying holons. It has 
been pointed out [15,16] that the physics of underdoped cuprates is captured 
by an effective U(l) gauge theory of massless fermions and charged scalar 
fields 

N 

C F = Vvw (dp - ia p) + |(<9 M - ia^) 6| 2 + V(\b\ 2 ). (41) 

(J—1 


Here b = ( 61 , 62 ) is a doublet of scalar fields representing the holons [15]. 

= 1 and generally ^ v a ,2 as a result of the velocity anisotropy; 
however, for simplicity we can let = v a ^ = 1- Since the spin and charge 
degrees of freedom are assumed to be separated, there is no Yukawa-type cou¬ 
pling between the fermion held and the scalar held. In the superconducting 
state, the boson 6 acquires a nonzero vacuum expectation value, i.e., (6) ^ 0. 
This nonzero (6) spontaneously breaks gauge symmetry of the theory and 
the gauge boson acquires a finite mass £ via Anderson-Higgs mechanism. 

In the context of high temperature superconductors the U(l) gauge held is 
introduced as a Lagrangian multiplier to impose local no-double occupancy 
constraint. It has no kinetic term ~ F' 2 W and its dynamics is obtained by 
integrating out the matter helds. If we only include the one-loop diagram in 
the vacuum polarization, we get D^ l {q 2 ) = g 2 vr(g 2 ) + £. As we have showed 
previously, the effect of additional scalar doublet is to shift N in the gauge 
boson vacuum polarization n(q 2 ) to N + 1, i.e., 71 (q 2 ) — (N + l)/8 |g|. Then 
the propagator for the gauge boson is 


d,J:p - k) 


_8_ (r _ k)A 

(N + l)(\p - k\ +rj) \ ^ ( p-k) 2 J ' 


(42) 


From the corresponding DS equation in the non-local gauge obtained above, 
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we found a critical gauge boson mass at A/£ c „ = 100. For small £, CSB 
occurs; while for £ > £ c „, CSB is suppressed. 

We now would like to discuss the long-range behavior of the AF correla¬ 
tion function. The AF spin correlation is defined as 

1 r rl 3 1 r 

(S + S -)„ = --J j—^Tr[a„(k)G 0 (k+p )], (43) 

where Go(k) is the fermion propagator. If the fermions are massless, then 

Go(k) = -2A, ( 44 ) 

and we have 

(S+S-) „ = -M (45) 

At p —> 0, (S' + S' - )o —s► 0, and the AF correlation is heavily lost. This is not a 

surprising result since our starting point is the resonating valence bond (RVB) 
picture proposed by Anderson [21], which is just a liquid of spin singlets and 
hence it has only short range AF correlation. This is not a satisfying situation 
because a long-range Neel order was observed in experiments shortly after 
the discovery of cuprate superconductors. 

However, even if we starts from an RVB ground state, it is still possible to 
obtain the long-range AF correlation because of the strong correlation nature 
of the Mott insulators. The strong correlation is reflected in the fact that 
double occupancy on a single lattice is completely inhibited due to the strong 
Coulomb repulsive force. After this local constraint and quantum fluctua¬ 
tions are taken into account, a strong U(l) gauge field emerges in the effective 
theory. This gauge held has important effect on physical properties since it 
can cause fermion condensation and give the originally massless fermions a 
finite mass. The AF spin correlation is expected to be significantly enhanced 
once the fermions become massive. To show this actually happens, we will 
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calculate the spin correlation function in the CSB phase (see Appendix for 
details). Although the dynamically generated fermion mass depends on the 
3-momentum, here, for simplicity, we assume a constant mass m for the 
fermions. This approximation is valid because we only care about the low- 
energy property and £(p 2 ) is actually a constant at p —> 0. The propagator 
for the massive fermion is 


G(k) 


— (7 • k + im ) 
k 2 + m 2 ’ 


(46) 


which leads to 

<S + S->o = 


p 2 + 4 m 2 

. m H-—:— arcsin 

4t r \ 2 |p | 



(47) 


This spin correlation behaves like —m/27r as p 0 and we have long-range 
AF correlation when CSB takes place. Therefore, strong fluctuations around 
the RVB ground state enhances the long-range AF spin correlations. 

We should emphasize that calculations based on perturbative expansions 
can not be used to obtain the long-range AF order. It might be argued 
that including higher order diagram can enhance the AF spin correlation. 
However, this argument is not right. If we include the gauge held while 
keeping the fermions massless, then the spin correlation is [37] 

< s+s_> - = -l2V(kT) lp ' ln (7) < 48) 

which damps at low momentum p —> 0. Rantner and Wen [37] used to claim 
that long-range AF correlation can be obtained by reexponentiating the spin 
correlation function [38]. This scenario is based on their previous statement 
[14] that the U(l) gauge held can not generate a Unite mass for fermions and 
hence is a marginal perturbation. This result is derived by considering only 
the one-loop correction of gauge held to the fermion self-energy. However, 
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CSB is a nonperturbative phenomenon and whether the gauge held generates 
a finite mass for the massless fermions can only be settled by investigating 
the self-consistent DS equation for the fermion self-energy. If the equation 
(12) does not have the nonlinear term in the denominator of the kernel, it is 
a linear equation and can not develop any genuine nontrivial solution. From 
the point of view of bifurcation theory, a linear operator has no bifurcation 
points those are necessary for a phase transition to take place. Once the 
nonperturbative effect is taken into account, the strong gauge held generates 
a finite fermion mass which breaks the chiral symmetry and gives rise to long- 
range AF order (Ref. [39] discussed the correspondence of CSB to AF order in 
another way). Actually, the formation of long-range AF order spontaneously 
breaks the rational symmetry of the system and generates a gapless spin 
wave excitation which corressponds to the massless Goldstone boson. These 
are hard to understand if we only include the gauge fluctuations without 
breaking any symmetry. Furthermore, the strong interaction of the gauge 
held with massless fermions of flavor 2 will unavoidably generate a finite 
fermion mass. 

Now we would like to discuss the application of our result to the interplay 
of various ground states in high temperature cuprate superconductors. It is 
well-known that the gauge boson mass £ is proportional to the superfluid 
density p, thus we can use £ to describe the superconducting order. Other¬ 
wise, we use CSB to describe the long-range AF order. Based on the fact 
that the superfluid density is proportional to doping concentration, we ob¬ 
tain a clear picture of the evolution of different orders upon increasing the 
doping concentration. At zero and low doping the gauge boson mass is zero 
or very small, so CSB and hence the AF order is present. When the doping 
concentration is larger than a critical value S cv the gauge boson acquire a 
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mass that is large enough to suppress the CSB and the AF order. Note that 
superconductivity begins to appear as the ground state of cuprate supercon¬ 
ductors at 5 sc which is less than 8 CV . Therefore, for S sc < 5 < S cv there is 
a coexistence of the AF order and the SC order in the bulk materials. Due 
to this coexistence, the length scale for AF order to appear should be larger 
than the vortex scale, which is consistent with STM experiments of Hoffman 
et al. [36]. 

When the external magnetic field is stronger than H c2 , the superfluid is 
completely suppressed, and, correspondingly, the gauge boson become mass¬ 
less. Then CSB reappears in the bulk material and gives a mass to the mass¬ 
less fermions. This mass provides a finite gap for the low energy fermions 
to be excited, thus at low temperature no fermionic excitations can be ob¬ 
served [15]. This causes the breakdown of the Wiedemann-Franz (WF) law 
in the normal ground state of cuprate superconductors [15,40]. Thus, based 
on spin-charge separation and CSB, we give a unified description for both 
the behaviors of AF spin correlation and the transport properties from a 
held theoretical point of view. This is the most noticeable advantage of our 
scenario comparing with so many other scenarios [41-47] those also address 
the problem of local AF order in vortex cores. 

6 INSTANTON EFFECT ON CSB 

Confinement is one unresolved problem in modern particle physics. A sem¬ 
inal paper written by Polyakov [22] has shed some light on this problem 
by studying a three-dimentional compact pure U(l) gauge theory (compact 
pure QED 3 ). In general, one can define an Abelian gauge theory on a two- 
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dimensional lattice which has the following action 

S = 5W 5] C 1 - cos F i,<*p) > ( 49 ) 

i,a/3 

with the field strength 

F \,a/3 -^i,a ~t~ -^-i+a,p -' 4 i+/3,a • (50) 

Here, the pairs (i,a) are used to denote the links between lattices, with i 
the beginning of a link and a its direction. If the vector potential A i a is 
defined to be a real number on its whole region, i.e., — oo < A^ a < +oo, 
the continuum limit of this action is just that of the usual U(l) gauge as 
presented in (1). However, a highly nontrivial physical effect emerges if the 
vector potential Hj jQ , has angular properties and hence is defined on a circle 
as —7r < Ai ta < 7r. Due to the periodicity of its action, such a field theory is 
called compact QED. 

Polyakov firstly considered the pure compact QED 3 without coupling 
matter fields to the gauge field. He found that instantons appear in this 
model as topological solutions of the Euclidean gauge field equations and lead 
to permanent confinement of static charges which is reflected by the area law 
for the Wilson integral. Compact QED 3 has attracted intense investigations 
in the past twenty years, initially as a simpler model to study quark confine¬ 
ment. Recently, compact QED 3 with matter fields has been used to model 
the physics of many strongly correlated electron systems [19,48]. However, 
although it is widely accepted that confinement is present in pure compact 
QED 3 , there is no consensus on the fate of permanent confinement when 
matter fields are included [49,50]. Comparing with compact QED 3 of scalar 
fields [49], the situation for compact QED 3 of massless fermions is particu¬ 
larly complicated because of the possibility of dynamical mass generation for 
the fermions. 
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Since compact QED 3 is originally defined on lattices, Monte Carlo nu¬ 
merical simulations are expected to provide important information on CSB, 
but they suffers from the notorious fermion sign problem. In this paper, 
we would like to analyze the chiral behavior using the DS equation method. 
To do this, we map the compact QED 3 onto a continuum theory and intro¬ 
duce the ultraviolet cutoff A keeping track of its lattice origin. As shown 
by Polyakov, the gauge field acquires a finite mass due to Debye screening 
caused by the instantons. We can use the mass £ of gauge field to describe 
the instanton effect and investigate the relation between CSB and instanton 
effect by solving the DS equation that consists of a massive gauge boson 
propagator. 

We have studied the relation of gauge boson mass and CSB in the context 
of high temperature superconductors in the last section. However, the critical 
gauge boson mass £ c „ is very small in the presence of additional scalar fields, 
due to the shift from N to N + 1 by the scalar doublet, and hence CSB can 
exist only for a small region of £. But the critical value £ c „ in the present 
case (compact U(l) gauge field coupled only to massless fermions) is much 
larger and there is a wide region of £ for CSB to take place. In the previous 
papers [5,6] addressing CSB in QED 3 the ultraviolet cutoff is provided by 
a = N/8 , which is kept fixed when the fermions flavor N is taken to infinity, 
because for momentum p > a the self-energy function damps rapidly. From 
'f| = £ 2 /« we know that the critical gauge boson mass is about £ c = a/2. The 
instanton effect can coexist with CSB for £ < a/2. If we couple a fermion 
of one flavor to compact gauge field, then CSB can coexist with instanton 
effect in a much wider region of £. 

The above result can be used to investigate the possible confinement to 
deconfinement transition in compact QED 3 because whether the fermions 
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have a finite mass is expected to affect the fate of permanent confinement 
[51]. Such a transition is no doubt of great importance in both particle 
physics and condensed matter physics, but beyond the scope of this paper. 

7 Summary and Discussion 

In this paper, we have discussed the effect of a finite mass of U(l) gauge boson 
on CSB and its physical implications. The gauge boson mass £ is reflected 
in the modification of the gauge field propagator, which appears in the DS 
equation of the fermion self-energy. The DS equation is nonlinear and hence 
hard to be solved. Iteration procedure is the most frequently used numer¬ 
ical calculation method, but it is not clear whether the iteration procedure 
leads to a convergent result or not. To avoid the problem brought by the 
convergency of iteration, we make use of bifurcation theory and parameter 
imbedding method to numerically investigate the DS equation. Adopting the 
Landau gauge and neglecting the wave function renormalization, we found 
a critical value £ c „ for the gauge boson mass that separates the CSB phase, 
for £ < £ c „, and chiral symmetric phase, for £ > £ c „. We then showed that 
including higher order corrections of the wave function renormalization does 
not change qualitatively the critical value £ c ,,. 

We then use our result to two physical systems, high temperature cuprate 
superconductors and compact QED 3 . If the gauge boson mass is generated 
via Anderson-Higgs mechanism in the superconducting state, the combina¬ 
tion of spin-charge separation and CSB provides a field theoretical description 
of the competition between the AF order and the SC order. As a compromise 
of this competition, there is a microscopic coexistence of these two orders in 
the bulk materials, which plays an essential role in explaining the local AF 
ordre in vortex states observed in neutron scattering and STM experiments. 
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When the periodicity of the gauge field is taken into account, the gauge boson 
acquires a finite mass via the instanton effect. Since whether the permanent 
confinement still exists in the presence of fermions depends on the fermion 
mass, our result can help us to investigate the confinement to deconfinement 
phase transition, which will be the subject of future study. 
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Appendix 


In this Appendix we give the details of calculating the spin correlation us¬ 
ing the propagator of the massive fermions. When CSB occurs the fermion 
propagator is 


G(k) = 


— (7 • k + im) 
k 2 + m 2 


(51) 


Then 


(S+S-) „= 4 / 


—l 


d 3 k 
(2t r) 3 
d 3 k 

W) 

r 

dt 


Tr [G(k)G(k + p)\ 

7 • k + im 7 • (k + p) + im 
k 2 + m 2 (k + p) 2 + m 2 

k ■ {k + p) — m 2 


-Tr 


d 3 k 


(27t) 3 [[k 2 + m?) (1 — t) + (( k + p) 2 + m 2 ) t] 
where we used the Feynman parameterization formula 

If 1 , 1 

— = / dt - -77. 

ab Jo [at + 6(1 — t) 2 } 


m 


(53) 
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After replacing k by k — pt and making Wick rotation, we have 


<S+S->„= Cdt < Cpk 


(54) 


Jo J (27t) 3 [k 2 + m? + p 2 t(l — t)} 2 

Using the properties of the T-function, we can integrate upon the momentum 
k and get 

1 


(m^ >r( i: r|{) f 1 * 


(4tt 


Since 


T(2) -to [m, 2 + p 2 t(l — t)] 2 

EH) 

r(2) 

r(i 


= —27T 2 


r (2) 


= 


we get 


( S + S ) 0 =- [ dt m 2 + p 2 t( 1 — t) 

2tc Jo L 


1 

47T 


p 2 + Am 2 
m H-—:— arcsm 

2 Ip 


p- 


p 2 + 4 m 2 


(55) 

(56) 

(57) 


(58) 


This is the AF spin correlation in the CSB phase. 
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Figure 1: The dependence of the critical number N c on 
log 10 (A/?/) in the Landau gauge. 
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Figure 2: The dependence of the critical number N c on log 10 (A/^) 
in the non-local gauge. 
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